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1.  Introduction 


The  experimental  determination  of  an  energetic  material’s  damage  potential  is  characterized  by 
the  detonation  wave  and  air  shock  wave  properties  produced  as  a  result  of  its  rapid 
decomposition.  The  key  detonation  wave  properties  needed  for  characterization  include: 
detonation  or  Chapman-Jouguet  (CJ)  pressure  Pcj,  detonation  wave  velocity  D,  reaction 
products  particle  velocity  ucj,  and  reaction  products  density  pCj.  The  ultra-high- speed  nature  of 
the  detonation  event  involves  nanosecond  timescale  chemical  reactions,  producing  temperatures 
between  2000-5000  K,  pressures  upwards  of  50  GPa,  and  detonation  wave  velocities  ranging 
from  4-10  mm/ps  (1-3).  Experimentally  measuring  these  properties  presents  a  difficult 
challenge  due  to  the  severe  nature  of  the  environment  and  the  temporal  scale  on  which  they  occur. 
Presently,  to  determine  the  desired  detonation  wave  properties,  a  multitude  of  experiments  (many 
conducted  using  kg-range  charges)  and  instrumentation  are  necessary.  These  can  include 
piezoelectric-pin  instrumented  detonation  velocity  experiments,  plate  dent  or  aquarium  test 
experiments,  and  photon  Doppler  velocimetry  experiments,  to  name  a  few.  Thus,  the  current 
detonation  characterization  protocol  ultimately  limits  the  development  of  novel  energetic 
formulations  as  it  tends  to  become  cost  prohibitive. 

Early  detonation  characterization  research  on  the  detonation  wave  reaction-zone  length  and 
pressure  was  first  suggested  by  Goranson  in  1946  (4).  He  determined  that  the  reaction-zone 
length  and  CJ  pressure  for  an  energetic  material  could  be  determined  through  free-surface 
velocity  measurements  of  metal  “flyer  plates”  positioned  adjacent  to  an  energetic  charge.  Duff 
and  Houston  later  tested  this  theory  by  characterizing  Composition  B  using  aluminum  flyer  plates 
of  various  thicknesses  (5).  Further  detonation  characterization  was  conducted  by  Holton,  who 
extrapolated  CJ  pressures  using  the  “aquarium  test,”  an  extension  to  Goranson’s  flyer-plate 
experiment,  in  which  the  flyer  plate  is  replaced  by  an  optically  transparent  material  (most 
commonly,  water)  (<5).  In  this  test,  a  streak  image  is  taken  parallel  to  the  axis  of  a  sympathetically 
detonated  cylindrical  explosive  charge  submerged  in  a  container  filled  with  distilled  water.  The 
recorded  streak  record  is  subsequently  analyzed  to  determine  the  time-dependent  position  of  the 
outward-propagating  shock  wave.  An  empirical  shock  wave  distance-versus-time  relationship  is 
applied  to  the  data  and  from  this,  the  shock  wave  velocity  at  the  water-charge  interface  may  be 
determined  through  its  differentiation  and  extrapolation  to  the  water-charge  interface.  From  this, 
a  detonation  pressure  is  calculated  using  Goranson’s  equation  with  the  known  Hugoniot  for  water 
and  known  energetic  material  and  water  densities  (7). 
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Further  aquarium  test  research  was  conducted  by  Rigdon  and  Akst  to  determine  the  optimum 
procedure  for  extrapolating  the  shock  wave  velocity  at  the  material  interface,  as  aquarium  test 
detonation  pressures  were  consistently  lower  than  published  values  (8).  A  host  of  empirical 
curvefitting  techniques  were  applied  to  a  range  of  data  sets  to  determine  the  optimum  treatment  of 
the  shock  wave  radius-versus-time  data.  They  concluded  that  a  linear  regression  fit  to  the  initial  2 
mm  of  shock  wave  travel  (7.2  cm  charge  diameter)  yielded  the  highest-accuracy  results  when 
compared  to  published  values. 

Conducting  the  air  shock  wave  characterization  of  energetic  materials  presents  a  similarly 
challenging  environment  due  to  the  pressures  and  temperatures  being  transferred  into  the 
surrounding  medium.  Air  shock  wave  characterization  properties  of  interest  include  the  shock 
wave  pressure  P,  shock  wave  velocity  U,  particle  velocity  u,  explosive  impulse  I,  and  more 
recently  of  interest,  fireball  temperatures  T  (9).  Air  shock  wave  characterization  is  traditionally 
conducted  using  kg-range  charges  instrumented  with  piezoelectric  or  piezoresistive  pressure 
transducers  positioned  at  specified  radial  standoff  distances.  These  transducers  measure  the 
shock  wave  time  of  arrival,  peak  shock  wave  pressure,  and  explosive  impulse,  and  allow  a  shock 
wave  velocity  to  be  determined  given  multiple  transducers.  Due  to  the  harsh  nature  of  the 
environment,  however,  measurements  are  rarely  taken  at  radial  standoffs  of  <5  charge  diameters, 
providing  little  data  suitable  towards  deriving  the  detonation  wave  or  air  shock  wave  properties  at 
the  energetic  material-air  interface. 

In  response  to  these  aforementioned  shortcomings,  a  laboratory-scale-based  modification  of  the 
aquarium  test  has  been  developed  to  allow  both  the  detonation  wave  and  the  air  shock  wave 
properties  at  the  energetic  material-air  interface  to  be  determined  from  a  single  experimental 
measurement.  The  optically  based  method  incorporates  high-speed  photography  and 
laboratory-scale  ( i.e .,  gram-range)  spherical  explosive  charges.  The  detonation  wave  and  air 
shock  wave  properties  for  an  energetic  material  of  interest  may  be  determined  from  material 
Hugoniots,  conservation  laws,  two  empirically  established  relationships,  and  a  single  film-based 
streak  record.  The  laboratory-scale -based  method  provides  an  accurate,  cost-effective, 
characterization  procedure  for  novel  energetic  materials  as  an  alternative  to  the  full-scale 
characterization  experiments  currently  used. 
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2.  Theory/Methodology 


A  shock  wave  is  fully  defined  through  knowledge  of  five  state  variables:  shock  pressure  P,  shock 
velocity  U,  particle  velocity  u,  density  p  or  specific  volume  v,  and  specific  internal  energy  e. 
These  variables  are  defined  by  the  conservation  of  mass  (equation  1),  momentum  (equation  2), 
and  energy  (equation  3)  jump  equations  across  the  shock  wave  (subscripts  0  and  1  refer  to  the 
states  before  and  after  the  shock’s  passage).  In  order  to  solve  for  the  shock  wave’s  state,  two 
additional  relationships  are  necessary.  The  first  relationship  is  an  equation  of  state  for  the 
material  of  interest.  The  desired  equation  of  state  would  relate  the  specific  internal  energy  as  a 
function  of  pressure  and  specific  volume,  i.e.,  e  =  f(P,  v).  By  combining  it  with  the 
conservation  of  energy  jump  equation  (equation  3),  the  specific  internal  energy  is  eliminated, 
leaving  the  well-known  Hugoniot  equation,  equation  4.  The  Hugoniot  equation  describes  the  loci 
of  pressure  states  attainable  after  a  shock  wave’s  passage  through  a  material  as  a  function  of  the 
material’s  specific  volume.  This  reduces  to  four  variables  (P,  U,  u,  and  v)  in  three  equations 
(equations  1,  2,  and  4).  Thus,  if  an  additional  relationship  may  be  determined  between  the  four 
remaining  variables,  the  system  is  closed. 


pi  _  U  -  u0  _  Vo 
p0  U  -  ui  v1 


Pi  -  P0  =  po(wi  -  U0)(U  -  Mo)  (2) 


ei  —  eo  — 


P\  U\  -  PqUq 
Po(U  -  Wo) 


(3) 


P  =  f(v)  (4) 

Material  Hugoniots  are  experimentally  measured  by  subjecting  a  material  to  a  series  of  known 
shock  compression  experiments.  A  total  of  six  Hugoniot  “planes,”  or  pairs,  may  be  determined 
based  on  the  four  remaining  state  variables:  P  —  U,  P  —  u,  P  —  v,  U  —  u,  U  —  v,  and  u  —  v.  Of 
these,  three  planes  have  been  found  to  be  of  great  importance  when  analyzing  shock  wave 
interactions:  U  —  u,  P  —  v,  and  P  —  a.  Thus,  adding  the  U  —  u  or  P  —  u  Hugoniot  relationship 
allows  the  shock  wave’s  state  to  be  fully  defined.  By  extending  these  equations  to  the  interaction 
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of  an  energetic  material  detonating  in  a  surrounding  fluid  medium,  each  material’s  “shocked” 
state  may  therefore  be  determined. 

Using  these  equations,  it  is  proposed  to  calculate  the  detonation  wave  and  transmitted  shock  wave 
properties  at  the  interface  of  an  energetic  material  detonating  in  air  through  the  sole  measurement 
of  the  air  shock  wave  velocity  U.  From  this  measurement,  the  following  seven  properties  are  to 
be  solved:  detonation  wave  velocity  D,  detonation  wave  pressure  Pcj ,  reaction  products  particle 
velocity  uCj,  reaction  products  density  pCj,  air  shock  wave  pressure  P,  air  particle  velocity  u, 
and  air  density  p  .  Thus,  a  system  of  seven  equations  is  necessary  to  define  the  detonation 
interaction. 

The  first  two  equations  implemented  are  the  conservation  of  mass  and  momentum  equations 
across  a  detonation  wave.  Here,  the  equations  are  written  for  the  case  of  a  one-dimensional 
planar  detonation  wave.  Conservation  of  mass  (equation  1)  is  rewritten  for  the  detonating 
energetic  material,  yielding  equation  5.  For  this,  the  density  ratio  is  the  reaction  products  density 
pcj  divided  by  the  energetic  material  pressing  density  pe,  detonation  velocity  D  is  substituted  for 
the  shock  wave  velocity  U,  and  the  reaction  products  particle  velocity  ucj  replaces  the  particle 
velocity  behind  the  detonation  wave  u  \ .  The  particle  velocity  ahead  of  the  detonation  wave  u0  is 
zero  as  the  material  is  yet  to  be  affected  by  the  detonation. 


Pcj  =  D 
pe  D  -  uCj 


(5) 


Detonation  wave  property  variables  are  additionally  substituted  into  the  conservation  of 
momentum  equation  (equation  2).  The  previous  conservation  of  mass  substitutions  remain  valid 
and  the  pressure  ahead  of  the  detonation  wave  is  assumed  to  be  zero  as  P0  Pcj .  Thus, 
conservation  of  momentum  across  the  detonation  wave  is  represented  by 


Pcj  =  peUcjD  .  (6) 

When  the  outward  propagating  detonation  wave  reaches  the  energetic  material-air  interface,  a 
transmitted  shock  wave  is  produced  and  driven  radially  outward  through  the  air  by  the  expanding 
energetic  material  reaction  products  (10).  As  a  result  of  this  interaction,  a  rarefaction  wave  is 
driven  rearward  through  the  reaction  products  toward  the  point  of  initiation  (2).  To  satisfy  the 
laws  of  conservation,  the  pressure  and  particle  velocity  at  the  energetic  material-air  interface 
must  be  equal  for  the  forward  propagating  transmitted  shock  wave  and  the  rearward  propagating 
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rarefaction.  Thus,  if  a  pressure-particle  velocity  (P  —  u)  Hugoniot  is  known  for  both  the 
energetic  material  reaction  products  and  air,  a  solution  for  the  pressure  and  particle  velocity  at  the 
interface  is  determinable. 

Presently,  there  exists  a  limited  collection  of  Hugoniot  data  for  energetic  material  reaction 
products  due  to  the  numerous  existing  explosives  and  their  multiple  densities  incorporated 
throughout  the  field.  In  response  to  this,  Cooper  developed  a  non-dimensional  empirical 
relationship  from  published  P  —  u  Hugoniot  data  for  energetic  material  reaction  products  (2).  He 
discovered  that  if  the  compilation  of  existing  Hugoniot  data  are  plotted  in  a  reduced  form,  i.  e. , 
reduced  pressure  P/Pcj  and  reduced  particle  velocity  u/ucj,  they  collapse  to  a  single  profile, 
represented  by  two  empirical  relationships  (2).  This  profile  is  divided  into  two  equations  based 
upon  the  reduced  pressure  value.  For  reduced  pressures  greater  than  0.08,  the  reduced  reaction 
products  Hugoniot  is  represented  by  equation  7  and  has  a  coefficient  of  determination  of  0.987. 

—  =  2.412-  1.7315  f— 1  +0.3195  f— 1  (7) 

Pcj  \UcjJ  \UcjJ 

For  reduced  pressures  less  than  0.08,  the  reduced  reaction  products  Hugoniot  is  represented  by 
equation  8  and  has  a  coefficient  of  determination  of  0.898. 


P 

P^i 


(8) 


For  the  present  case,  the  reduced  pressure  is  anticipated  to  be  extremely  low  due  to  the  large 
shock  impedance  discrepancy  between  the  energetic  material  and  air.  Shock  impedance  Z  is 
defined  as  the  product  of  a  material’s  density  and  shock  wave  velocity  (2),  equation  9.  The 
energetic  material  impedance  Ze  is  defined  as  the  product  of  its  pressing  density  pe  and  its 
detonation  velocity  D  (equation  10),  whereas  air’s  shock  impedance  Zair  is  defined  as  the  product 
of  its  ambient  density  pair  and  the  transmitted  shock  wave  velocity  U  (equation  11).  The 
energetic  material  pressing  density  is  ~1500  times  greater  than  the  density  of  air.  If  it  is  assumed 
that  the  detonation  wave  velocity  and  air  shock  wave  velocity  are  of  approximately  the  same  order 
of  magnitude,  the  detonation  wave  interaction  with  air  results  in  a  transmitted  shock  wave  having 
a  pressure  less  than  the  detonation  pressure  Pcj  and  a  particle  velocity  greater  than  ucj 
according  to  the  material  Hugoniots  (figure  1).  Knowledge  of  the  air  shock  wave-particle 
velocities  relationship  is  required  to  determine  how  low  this  reduced  pressure  will  be.  It  should 
be  noted  that  each  material’s  impedance  is  not  constant  and  that  both  will  increase  slightly  with 
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increasing  pressure.  However,  for  practical  purposes,  the  shock  impedances  may  be  considered 
constant  across  a  reasonable  range  of  interest  (2). 


Z  =  pU  (9) 

=  peD  (10) 


Zair  PairU 


(11) 


Particle  velocity,  u 

Figure  1.  Interface  interaction  on  the  P  —  u  Hugoniot  plane. 


Previous  research  has  shown  that,  for  most  materials,  a  linear  relationship  exists  between  the 
shock  wave  and  particle  velocities,  or  the  U  —  u  Hugoniot  (figure  2).  The  linear  relationship  is 
composed  of  two  empirically  measured  coefficients:  Cq,  defined  as  the  bulk  sound  speed  (though 
it  has  no  physical  meaning)  and  s,  corresponding  to  the  U  —  u  Hugoniot  slope,  equation  12.  The 
U  —  u  Hugoniot  coefficients  for  air  are  well  established  from  previous  research  by  Deal  who 
conducted  plate  impact  experiments  using  explosively  driven  flyer  plates.  In  his  research,  he 
determined  air’s  Hugoniot  coefficient  Cq  to  be  0.2375  km/s  and  s  to  be  1.0575  for  measured 
shock  wave  velocities  up  to  4.5  km/s.  A  higher  shock  wave  velocity  is  expected,  however,  air’s 
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Hugoniot  is  assumed  to  remain  linear  throughout  the  measurement  regime.  Given  that  the  air 
shock  wave  velocity  is  linearly  related  to  the  particle  velocity  (equation  12)  and  the  energetic 
material  pressing  density  is  ~1500  times  greater  than  the  density  of  air,  an  extremely  large 
impedance  mismatch  between  the  materials  results,  and  ultimately,  a  much  lower  pressure  is 
produced  (2).  Thus,  the  reduced  pressure  may  be  anticipated  to  fall  below  the  0.08  established 
threshold,  and  equation  8  is  used  for  solving  the  system.  Lastly,  it  should  be  noted  that  the  U  —  u 
Hugoniot  for  air  is  used  in  place  of  the  P  —  p  Hugoniot  as  it  has  previously  been  well 
characterized. 


Particle  velocity,  u 
Figure  2.  U  —  u  Hugoniot  plane. 


U  =  Co  +  su  (12) 

As  shown  by  figure  1,  an  equation  representing  the  P  —  u  Hugoniot  plane  for  air  is  required  to 
determine  the  shock  wave  pressure  and  particle  velocity  solution  to  the  energetic  material-air 
interaction.  This  is  obtained  by  combining  the  conservation  of  momentum  equation  (equation  2) 
and  air’s  U  —  u  Hugoniot  relation,  yielding  equation  13.  As  with  conservation  of  momentum 
across  the  energetic  material  detonation  wave  (equation  6),  the  shock  wave  velocity,  shock  wave 
pressure,  and  particle  velocity  ahead  of  the  shock  wave  (U0,  P0,  and  u0  in  equation  2, 
respectively)  are  zero,  as  the  air  is  yet  to  be  affected  by  the  detonation  and  P0  -C  P.  Thus,  given 
equations  8  and  13,  a  solution  for  the  shock  wave  pressure  and  particle  velocity  at  the  energetic 
material-air  interface  is  determinable  by  plotting  the  two  P  —  u  Hugoniots. 
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(13) 


P  Pair  P{}  U  ~t~  Pair 

The  final  equation  required  to  define  the  air  shock  wave  is  the  conservation  of  mass  written  for 
across  a  shock  wave,  equation  14.  For  this,  the  density  ratio  is  the  air  shock  wave  density  p 
divided  by  the  atmospheric  density  pair,  the  shock  wave  particle  velocity  u  is  the  particle  velocity 
behind  the  detonation  wave  U\,  and  the  shock  wave  velocity  U  remains.  As  before,  the  shock 
wave  velocity,  shock  wave  pressure,  and  particle  velocity  ahead  of  the  shock  wave  (U0,  Po,  and 
Uq,  respectively)  are  zero,  as  the  air  is  yet  to  be  affected  by  the  detonation  and  P0  <C  P. 


Pair  U  U 

The  final  equation  required  to  close  the  system  is  an  empirical  relationship  developed  by 
Cooper  (2)  that  relates  the  reaction  products  density  to  the  energetic  material  pressing  density, 
equation  15.  This  is  chosen  in  place  of  a  P  —  p  Hugoniot  for  the  reaction  products  as  no 
all-energetic-material-encompassing  Hugoniot  exists.  Equation  15  is  shown  to  be  valid  for 
densities  ranging  between  2e-3-  9  g/cm3. 


pCj  =  1.386 p°e96  (15) 

By  solving  the  system  of  equations,  the  desired  detonation  wave  properties  ( D ,  Pcj,  ucj,  and 
pcj )  and  air  shock  wave  properties  at  the  interface  (P,  u,  and  p)  are  determinable  for  an  energetic 
material  of  interest  by  measuring  the  air  shock  wave  velocity  U  at  the  energetic  material-air 
interface.  The  system  solution  is  represented  by  equations  16-22: 


U-C0 

u  = - 

s 


(16) 


P 


PairC qU  T  Pair^^~ 


(17) 


pCj  =  1.386pg'96 


(18) 


8 


In 


uCj  =  exp 


pu8.71 

235  PcJP* 

\  PCJ  -  Pe  J 
10.71 


(19) 


D  = 


ucjPcj 
PCJ  —  Pe 


(20) 


PCJ  =  PeUcjD  . 


(21) 


p  = 


UPa 


U -u 


(22) 


3.  Experimental  Methods 


3.1  Spherical  Charges 

Centrally  initiated  laboratory-scale  spherical  charges  are  used  to  eliminate  the  critical  diameter 
effect  associated  with  energetic-material  detonations  (11).  As  previously  mentioned,  the 
conservation  equations  developed  are  written  for  the  case  of  a  planar  shock  wave  interaction. 
However,  previous  research  has  demonstrated  that  the  radial  divergence  effects  associated  with 
spherical  detonation  waves  are  insignificant  for  shock  wave  curvatures  that  are 
orders-of-magnitude  larger  than  the  detonation-reaction-zone  thickness  (12-15).  Thus,  the 
planar  equations  remain  valid  for  the  present  explosive  charge  configuration. 

Spherical  charges  are  composed  of  matching  hemispheres  having  a  nominal  equivalent  mass  of 
0.90  g.  Each  hemisphere  possesses  a  centrally  positioned  half  right-circular-cylinder  void  to 
accommodate  a  RP-3  exploding  bridgewire  (EBW)  detonator  (29  mg  pentaerythritol  tetranitrate 
[PETN]),  when  assembled  as  shown  in  figure  3.  Charges  are  composed  of  Class  V  RDX  pressed 
to  a  density  of  1.77  or  1.63  g/cm3.  The  assembled  charge  masses  m  and  their  respective  densities 
pe  are  tabulated  in  table  1. 
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Figure  3.  Fully  assembled  spherical  charge  with 
RP-3  EBW  detonator. 


Table  1.  Experimental 

charges  detonated. 


Charge 

Pe 

(g/cm3) 

m 

(g) 

11318-3 

1.77 

1.809 

11318-4 

1.77 

1.814 

12040-2 

1.63 

1.779 

12040-4 

1.63 

1.749 

12052-1 

1.63 

1.726 

12052-3 

1.63 

1.758 

3.2  Charge  Imaging 

A  Beckman  &  Whitley,  Model  770  synchronized  film-based  streak  camera  is  used  to  record  the 
explosively  driven  air  shock  wave  radial  expansion.  The  streak  camera  has  a  slit  width  of  0. 1  mm 
and  a  write  speed  of  8  mm//is.  A  90  cm  // 6.3  field  lens  is  used  to  relay  a  focused  image  onto  the 
streak  camera  film  plane.  Figure  4  shows  the  experimental  setup  used  for  the  present 
experiments.  The  explosive  fireball  and  ionization  of  the  surrounding  air  due  to  the  shock  wave 
provide  self-illumination  of  the  record  from  which  radial  expansion  rates  are  measurable.  A 
geometric  correction  factor  is  not  implemented  into  the  present  measurements  as  the  angles 
associated  with  the  radial  expansion  measurements  are  extremely  small  (i.e.,  O(mrads)) , 
therefore  allowing  the  small-angle  approximation. 
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Figure  4.  Experimental  imaging  setup. 


4.  Results  and  Discussion 


4.1  Experimental  Results 

Prior  to  conducting  an  experiment,  a  static  image  of  a  calibration  object  (here,  a  dimensional 
scale)  was  recorded  to  determine  image  magnification,  and  ultimately,  a  calibration  length  scale. 
Spherical  charges  were  oriented  with  the  cylindrical  detonator  axis  collinear  to  the  slit  to  ensure 
the  transferred  air  shock  wave  measured  was  produced  by  a  steady-state  detonation  and  normal  to 
the  slit  axis  (figure  5).  As  shown,  two  fiducial  wires  were  oriented  perpendicular  to  the  slit  axis 
at  the  streak  camera  image  plane  to  allow  charge  orientation  and  position  to  be  determined 
between  static  and  streak  images.  The  explosive  fireball  illumination  and  air  shock  wave 
ionization  were  used  for  film  exposure.  This  produced  an  acceptable  streak  record  (50  /xs  total), 
as  only  the  initial  130  ns  were  necessary  for  the  present  research  (figure  6). 

Streak  records  were  analyzed  on  a  vertical-beam  optical  comparator  (S-T  Industries,  Inc.  Model 
4600)  with  20X  objective  lens.  Measurements  were  taken  at  0.05-mm  intervals  along  the 
temporal  axis,  corresponding  to  a  radial  shock  wave  measurement  taken  every  6.25  ns.  Data 
were  subsequently  transformed  into  dimensional  shock  wave  radius  and  time  coordinates  using 
the  calibration  length  scale  and  known  write  speed.  A  best-fit  linear  correlation  was  determined 
for  each  data  set  according  to  equation  23,  which  consists  of  the  shock  wave  radius  R,  air  shock 
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Figure  5.  Static  image  of  charge-slit 
orientation. 


Figure  6.  Sample  streak  record  at  8  mni///s 
writing  speed  (time=ordinate  and 
radius=abscissa)  . 


wave  velocity  U,  time  t,  and  the  spherical  charge  radius  r0.  Only  the  initial  1.2  mm  of  radial 
shock  wave  travel  were  measured  according  to  the  results  of  Rigdon  and  Akst  ( 8 ).  From  their 
research,  it  was  determined  that  a  linear  correlation  fit  to  the  initial  2  mm  of  radial  shock  wave 
travel  yielded  the  highest  accuracy  for  determining  detonation  pressures  from  aquarium 
experiments.  However,  their  research  incorporated  energetic  charges  having  diameters  of  7.2  cm. 
Using  the  explosive  scaling  laws,  it  was  determined  that  only  the  first  1.2  mm  of  radial  shock 
wave  travel  were  necessary  for  accurate  reduction  of  the  current  data.  Resulting  data  and  the 
average  linear  regression  fit  are  shown  in  figures  7  and  8. 


R  =  Ut  +  r0 


(23) 
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Time,  t  [f^s] 


Figure  7.  Shock  wave  radius- versus-time  data  for  RDX  at  1.77 
g/cm3  with  average  linear  regression  tit. 


Figure  8.  Shock  wave  radius- versus-time  data  for  RDX  at  1.63 
g/cm3  with  average  linear  regression  fit. 
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The  shock  wave  velocity  was  determined  from  the  best-fit  linear  correlation  to  the  data  and  input 
into  a  MATLAB  script  written  to  calculate  the  desired  characterization  properties.  The  script 
allows  the  user  to  input  the  energetic  material  pressing  density  pe  and  the  ambient  density  of  air 
pair  at  the  time  of  experiment.  Detonation  wave  and  air  shock  wave  properties  were  subsequently 
solved  for  each  experiment  according  to  equations  16-22.  Experimental  results  are  tabulated  in 
tables  2  and  3. 


Table  2.  Calculated  detonation  wave  and  air  shock  wave  properties  for  RDX  at  1.77  g/cm3. 


Explosive 

Pair 

(g/cm3) 

U 

(mm//is) 

D 

(nim//is) 

Pcj 

(GPa) 

PCJ 

(g/cm3) 

UCJ 

(mm//is) 

P 

(MPa) 

P 

(g/cm3) 

u 

(mm//ts) 

Published 

- 

- 

8.69  (16) 

33.6  (17) 

2.34  (2) 

- 

- 

- 

- 

11318-3 

1.202e-3 

9.18 

8.73 

35.3 

2.40 

2.28 

93.3 

1.532e-2 

8.46 

11318-4 

1.202e"3 

8.96 

8.51 

33.6 

2.40 

2.23 

88.9 

1.517e-2 

8.25 

Average 

1.202e"3 

9.07 

8.62 

34.4 

2.40 

2.26 

91.1 

1.524e"2 

8.36 

%  Deviation 

- 

- 

0.81 

2.53 

2.56 

- 

- 

- 

- 

Table  3.  Calculated  detonation  wave  and  air  shock  wave  properties  for  RDX  at  1.63  g/cm3. 


Explosive 

Pair 

(g/cm3) 

U 

(mm//js) 

D 

(mm///s) 

Pcj 

(GPa) 

PCJ 

(g/cm3) 

UCJ 

(mm//j,s) 

P 

(MPa) 

P 

(g/cm3) 

u 

(mm//is) 

Published 

- 

- 

8.34  (16) 

28.3  (17) 

2.16(2) 

- 

- 

- 

- 

12040-2 

1.207e"3 

8.85 

8.40 

30.4 

2.22 

2.22 

87.0 

1.504e-2 

8.14 

12040-4 

1.207e-3 

9.06 

8.61 

31.9 

2.22 

2.28 

91.2 

1.519e-2 

8.34 

12052-1 

1.213e"3 

9.16 

8.71 

32.7 

2.22 

2.31 

93.7 

1.543e-2 

8.44 

12052-3 

1.206e-3 

8.74 

8.30 

29.7 

2.22 

2.19 

84.8 

1.506e-2 

8.04 

Average 

1.208e-3 

8.95 

8.51 

31.2 

2.22 

2.25 

89.2 

1.518e-2 

8.24 

St.  Dev. 

3.202e"6 

0.19 

0.19 

1.4 

0.00 

0.05 

4.0 

1.794e"4 

0.18 

%  Deviation 

- 

- 

1.98 

10.2 

2.78 

- 

- 

- 

- 

The  percent  deviation  between  average  and  published  values  was  determined  using  equation  24. 
As  shown  by  tables  2  and  3,  the  detonation  wave  properties  (D,  PCj,  and  pcj)  were  all  within  3% 
of  their  published  values  with  the  exception  of  the  detonation  pressure  for  the  1.63  g/cm3  charges, 
which  was  found  to  be  10.2%.  It  should  be  noted,  however,  that  the  published  detonation 
pressure  for  RDX  pressed  to  1.63  g/cm3  was  empirically  derived  by  Keshavarz  as  an 
experimentally  measured  value  was  unavailable  (77).  Thus,  the  presently  calculated  deviation 
would  be  affected  by  the  inherent  empirical  error  associated  with  the  published  value. 
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%  deviation 


|  average  -  published\ 
published 


(24) 


An  analysis  was  performed  to  determine  the  inherent  uncertainty  within  the  calculated  data  due  to 
the  experimental  uncertainty  present  and  empirical  relationships  employed  according  to 
Moffat  (18).  The  uncertainty  in  measuring  the  air  shock  wave  velocity  U  was  calculated  from  the 
experimental  uncertainties  introduced  as  a  result  of  measuring  the  initial  charge  radius  (Ar0  = 
0.005  mm),  shock  wave  radius  (A R  =  0.005  mm),  and  time  (At  =6.25e_5/is).  This  led  to  an 
uncertainty  of  0.09  mm///s  for  the  air  shock  wave  velocity  U  and  0.08  mm//is  for  the  air  particle 
velocity  u.  The  uncertainty  in  the  calculated  reaction  products  density  pcj  was  determined  to  be 
6.5c  3  g/cm3,  resulting  from  the  experimental  uncertainty  in  the  pressing  density  measurement. 
These  uncertainties,  in  addition  to  the  uncertainties  in  the  atmospheric  density  at  the  time  of  test, 
produced  an  uncertainty  of  2  MPa  for  the  calculated  air  shock  wave  pressure  P.  The  reaction 
products  particle  velocity  ucj  possessed  an  uncertainty  of  0.02  mm//js  and  the  detonation  wave 
velocity  D  had  an  uncertainty  of  0.14  mm/^s.  Lastly,  the  detonation  pressure  Pcj  uncertainty 
was  determined  to  be  0.025  GPa.  Thus,  for  all  properties  calculated,  the  experimental  uncertainty 
was  determined  to  be  less  than  2%  of  their  respective  values. 

4.2  Historical  Data 

The  proposed  theory  was  applied  to  previous  research  incorporating  variants  of  the  aquarium  test 
to  determine  the  detonation  wave  and  air  shock  wave  properties  of  an  energetic  material.  The 
research  allowed  the  present  method  to  be  applied  to  a  range  of  full-scale  (i.e.,  kg-range) 
energetic  charges  across  various  atmospheric  conditions  (10,  19,  20).  The  tests  consisted  of 
experimentally  measuring  multiple  variables  for  a  given  explosive  as  well  as  having  a  priori 
knowledge  of  the  energetic  material’s  reaction  products  Hugoniot  or  detonation  variables.  In 
each  case,  the  measured  air  shock  wave  velocity  at  the  interface  and  atmospheric  density  were 
input  into  the  MATLAB  script,  allowing  all  other  detonation  wave  and  air  shock  wave  properties 
to  be  calculated.  The  results  for  several  cases  are  reported  in  table  4.  Air  shock  wave  density  at 
the  interface  was  not  reported  in  any  references  and  was  therefore  disregarded  in  the  present 
results.  The  developed  theory  was  able  to  predict  the  published  values  within  6%  (equation  25) 
for  all  properties,  while  the  majority  were  predicted  within  3%,  as  shown  by  the  tabulated  results. 
Thus,  the  theory’s  accuracy  was  extended  using  historical  characterization  research.  This  also 
demonstrates  the  theory’s  utility  as  the  only  properties  required  to  perform  a  fundamental 
detonation  characterization  on  an  energetic  material  of  interest  are  the  pressing  density, 
atmospheric  density,  and  the  initial  shock  wave  velocity. 
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%  deviation 


\calculated  -  published\ 
published 


*  100 


(25) 


Table  4.  Calculated  detonation  wave  and  air  shock  wave  properties  from  previous  experimental  research. 


Explosive 

Pe 

(g/cm3) 

Pair 

(g/cm3) 

U 

(mm//is) 

D 

(mm//is) 

Pcj 

(GPa) 

PCJ 

(g/cm3) 

P 

(MPa) 

u 

(mm//is) 

HMX/lnert  (79) 

(95/5) 

1.781 

1.068e-2 

7.42 

8.73 

33.5 

2.34 

53.0 

6.53 

Calculated 

- 

- 

- 

8.60 

34.5 

2.41 

53.8 

6.79 

%  Deviation 

- 

- 

- 

1.5 

3.0 

2.9 

1.5 

4.0 

HMX/T’JM  T/Inert  (19) 

(68/30/2) 

1.776 

3.785e-3 

7.84 

8.21 

31.2 

2.40 

21.2 

7.13 

Calculated 

- 

- 

- 

8.26 

31.7 

2.41 

21.3 

7.18 

%  Deviation 

- 

- 

- 

0.6 

1.7 

0.4 

0.5 

0.7 

TNT  (10) 

1.636 

0.97  le-3 

7.34 

6.93 

20.7 

2.15 

47.8 

6.71 

Calculated 

- 

- 

- 

6.80 

20.0 

2.22 

47.8 

6.71 

%  Deviation 

- 

- 

- 

1.9 

3.4 

3.2 

0 

0 

Composition  B  (10) 

(Grade  A) 

1.717 

0.947e"3 

8.67 

7.98 

29.5 

2.34 

65.5 

7.97 

Calculated 

- 

- 

- 

8.05 

29.2 

2.33 

65.5 

7.98 

%  Deviation 

- 

- 

- 

0.9 

1.0 

0.4 

0 

0.1 

Cyclotol  (10) 

(77/23) 

1.752 

0.936e-3 

8.64 

8.27 

31.3 

2.38 

64.2 

7.94 

Calculated 

- 

- 

- 

8.01 

29.5 

2.37 

64.2 

7.95 

%  Deviation 

- 

- 

- 

3.1 

5.8 

0.4 

0 

0.1 

Octol  (10) 

(77.6/22.4) 

1.821 

0.943e-3 

8.86 

8.49 

34.2 

2.46 

68.1 

8.14 

Calculated 

- 

- 

- 

8.23 

32.2 

2.46 

68.1 

8.15 

%  Deviation 

- 

- 

- 

3.1 

5.8 

0 

0 

0.1 

HBX-1  (20) 

1.712 

1.253e-3 

7.61 

7.31 

22.0 

2.18 

67.0 

6.97 

Calculated 

- 

- 

- 

7.22 

22.4 

2.21 

66.5 

6.97 

%  Deviation 

- 

- 

- 

1.2 

1.8 

1.4 

0.7 

0 
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5.  Conclusion 


A  laboratory-scale -based  modification  of  the  aquarium  test  has  been  provided  for  the  detonation 
wave  and  air  shock  wave  characterization  of  energetic  materials.  The  technique  measures  the 
radial  shock  wave  expansion  rate  produced  at  the  energetic  material-air  interface  from  the 
detonation  of  laboratory-scale  spherical  charges.  Air  shock  wave  velocities  were  measured  by 
fitting  a  regression  line  to  the  first  1 .2  mm  of  streak  record  data.  From  this  measurement,  a 
procedure  was  developed  to  determine  fundamental  detonation  wave  and  air  shock  wave 
properties  for  an  energetic  material  of  interest  through  use  of  conservation  laws,  material 
Hugoniots,  and  two  empirically  established  relationships.  The  theory  was  verified  using  pressed 
RDX  spheres  and  historical  aquarium  test  data. 

The  developed  theory  allows  a  fundamental  characterization  to  be  performed  from  a  single 
laboratory-scale  experiment  requiring  minimal  energetic  material,  offering  a  highly  advantageous, 
cost-saving  approach  for  the  screening  of  novel  energetic  formulations.  Thus,  it  provides  an 
alternative  to  the  multitude  of  full-scale  experimental  characterization  techniques  currently 
employed  during  the  developmental  phase  of  novel  energetic  formulations.  Presently,  to 
determine  the  desired  detonation  wave  (velocity  D,  pressure  PCj,  particle  velocity  uCj,  and 
density  pCj )  and  air  shock  wave  (pressure  P  and  particle  velocity  u)  properties,  a  multitude  of 
experiments  (many  conducted  using  kg-range  charges)  and  instrumentation  are  necessary. 
However,  as  shown  by  the  experimental  and  historical  data,  the  present  theory  only  requires  a 
single  experimental  measurement  to  obtain  all  of  the  above-named  properties.  The  theory’s 
utility  is  apparent  as  it  provides  fundamental  energetic  material  detonation  wave  and  air  shock 
wave  (at  the  energetic  material-air  interface)  characterization  properties  from  only  a  few  grams  of 
material.  Additional  testing  is  necessary,  however,  to  validate  the  theory  for  a  wider  range  of 
energetic  formulations.  Energetics  exhibiting  large  reaction-zone  thickness,  critical  diameter,  and 
energy  fluence  should  be  tested  to  verify  the  technique’s  validity.  Lastly,  it  is  recommended  that 
numerical  modeling  be  performed  to  ensure  that  the  assumed  second-order  spherical  divergence 
effects  are  indeed  negligible  for  the  proposed  laboratory-scale  charges. 
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List  of  Symbols,  Abbreviations,  and  Acronyms 


Go 

D 

EBW 

e 

HBX-1 

HMX 

/ 

m 

P 

Pcj 

PETN 

Po 

RDX 

s 

subscript 0 
subscript 
T 
To 

TNT 

U 

u 

UCJ 

V 

Pair 

Pcj 

Pe 

Ws 

Z 

7  , 

^  air 

ze 


air  bulk  sound  speed 
detonation  wave  velocity 
exploding  bridgewire 
specific  internal  energy 

pourable  mixture  of  RDX,  TNT,  and  aluminum 

cyclotetramethylene  tetranitramine 

explosive  impulse 

explosive  charge  mass 

shock  wave  pressure 

detonation  wave  pressure 

pentaerythritol  tetranitrate 

atmospheric  pressure 

cyclotrimethylene  trinitramine 

air  U  —  u  Hugoniot  slope 

properties  in  front  of  shock  wave 

properties  behind  shock  wave 

temperature 

atmospheric  temperature 
trinitrotoluene 
shock  wave  velocity 
particle  velocity 

reaction  products  particle  velocity 

specific  volume 

ambient  air  density 

reaction  products  density 

energetic  material  pressing  density 

streak  camera  write  speed 

shock  impedance 

air  shock  impedance 

energetic  material  shock  impedance 
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